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The  following  summary  is  a  progress  report  for  the  research  J 

currently  carried  out  under  AFOSR  Contract  #  AFOSR  81-0020.  The 
areas  covered  in  this  report  are  1) mathematical  theory  of  queer  — 
differential  equations  (QDE),  2)  universal  solutions  in  multi- 


£ 
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dimensional  diffusion  equations,  3)  exact  integrals  of  the  Emden-on- 
Fowler  equation,  4)  new  results  in  the  theory  of  turbulent  self- 


n/ 


diffusion  and  5)  mathematical  theory  of  the  essential  spectrum  in+y  Codes 

.Avail  and/or 


magnetohydrodynamics . 


1.  The  Theory  of  Queer  Differential  Equations  (QDE 1 s) 


Dist. 


tl 


Special 


Since  1974,  the  use  of  QDE's  in  the  theory  of  adiabatic  and 
diffusing  plasma  media  has  proved  very  useful.  More  recently, 


strong  interest  has  been  placed*  on  extending  the  use  of  QDE's  to 
other  problems  of  computation."  1  fluid  interest  where  there  are 
multiple  time  scales  such  that  they  can  be  separated  sufficiently 
to  allow  quasi-static  evolution  of  part  of  the  solution.  Large 
scale  meteorological  computations  are  envisioned  as  candidates 
for  such  an  application  of  QDE  techniques.  Numerical  algorithms 
have  been  developed  which  accelerate  large  diffusion  codes  and  the 
convergence  properties  of  the  numerical  schemes  depend,  to  a  large 
extent,  on  as  yet  incomplete  understanding  of  the  theory  of  QDE's. 
Recently,  research  in  the  mathematical  nature  of  QDE's  has  been 
pursued  by  H.  Grad,  P.  Laurence  and  E.  Stredulinsky .  We  present 
here  a  progress  report  of  this  work. 
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Equations  of  this  type  arise  in  the  work  of  H.  Grad  in  the  theory  of 


for 


ip  =  0  on  6ft,  ifi*(0)=0,  i|>M|ft|  =  1 


(2) 


In  other  words  ^(x)  is  zero  at  its  infimum,  and  one  at  its  supremum. 
Results  already  announced  in  the  last  interim  report  are  that  T(^) 
has  a  minimizer  in  a  function  class  that  incorporates  the  boundary 
conditions  (2). 

Several  regularity  results  have  been  established  for  this  pro¬ 
blem,  the  most  important  being  that  ip*  is  a  Lipschitz  continuous 
function  of  V  and  that  is  bounded  below  in  terms  of  an  expres- 

{*=t} 

sion  that  involves  only  the  measure  of  the  level  set  {^=t}  (see  [3]). 

It  rapidly  became  apparent  that  in  order  to  obtain  initial  re¬ 
sults  more  quickly,  at  least  initially,  it  would  be  necessary  to 
separate  questions  concerning  the  structure  of  the  level  sets  and  cri¬ 
tical  points  of  the  minimizer  from  its  differentiability*  (of  whatever 

*.  , 

order) .  For  this  reason,  an  auxiliary  variational  problem  is  intro¬ 
duced  [31  for  convex  ft  that  minimizes  T (xp )  among  functions  where 
level  sets  are  convex.  There  is  strong  evidence  that  a  minimizer  to 
this  auxiliary  problem  provides  a  solution  to  the  queer  differential 
equation  1) ,  even  though  the  admissibility  class  of  trial  functions 
has  been  narrowed.  All  details  regarding  this  are  not  complete  as  of 
this  writing  although  many  are  contained  in  [41.  In  any  case  this  is 
why  aside  from,  we  think,  being  interesting  in  their  own  right,  esti¬ 
mates  for  the  auxiliary  problem  will  carry  over  to  the  original  problem. 

For  the  auxiliary  problem  we  have  introduced  an  approximate  prob¬ 
lem  by  modeling  the  QDE  term  through  a  finite  difference.  The  nth 
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approxiroation  gives  rise  to  a  free  boundary  problem  with  n  free 
boundaries  which  is  a  topic  that  is  in  and  of  itself  on  the  fron¬ 
tiers  of  work  in  free  boundaries  (e.g.,  in  the  work  of  Luis 
Caffarelli  and  Avner  Friedman) .  It  can  be  shown  that  the  approxi¬ 
mate  problem  is  harmonic  in  the  regions  contained  between  the  free 
boundaries  and  that  certain  jump  conditions  on  the  gradient  of  the 
function  are  satisfied  across  the  free  boundaries.  A  combination  of 
the  use  of  these  jump  conditions  with  novel  results  on  V(t)  for 
harmonic  functions  permits  one  to  conclude  that  the  solution  to  our 
approximate  QDE  problem  is  superharmonic,  which  we  feel  is  itself 
a  significant  achievement.  Moreover,  as  the  property  of  being  super¬ 
harmonic  is  preserved  under  weak  limits  as  n  t  °°,  it  is  expected 
that  the  superharmonicity  result  will  carry  over  to  the  full  auxiliary 

problem  arid  thus  via  the  remarks  above  to  the  full  QDE. 

*  •»  *  * 

*  . . 

As  we  have  exhibited  a  strong 'interplay  between  understanding 
certain  free  boundary  problems  for  harmonic  functions  and. queer 
differential  equations,  an  interaction  between  us  and  L.  Caffarelli 
and  A.  Friedman  has  been  sparked  which  has  already  led  to  a  few 
new  theorems  which  will  also  be  contained  in  [4].  Furthermore, 
current  work  of  theirs  [5]  for  the  case  where  ft  is  not  convex  pro¬ 
vides  a  significant  boost  to  the  idea  of  attempting  a  similar  approxi¬ 
mation  procedure  by  an  n- she 11  free  boundary  problem  in  the  nonconvex 
case*  Queer  differential  equations  give  rise  to  many  new  and  in¬ 
teresting  questions  about  free  boundary  problems,  thus  we  think  they 
are  a  source  for  challenging  problems  for  years  to  come  in  this  area. 
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It  was  claimed  by  Harold  Grad  long  ago  that  "plasma  physics 
continually  gives  rise  to  new  and  interesting  mathematical  struc¬ 
tures."  The  relationship  we  have  sketched  between  queer  differen¬ 
tial  equations  and  free  boundary  problems  provides  another  example 
of  the  truth  of  Grad's  conjecture. 
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2.  Universal  Solutions  in  Multidimensional  Diffusion  Equations 
This  work  is  presently  being  carried  out  by  P.  Rosenau  and 
represents  a  natural  extension  of  his  work  in  non-linear  thermal 
and  diffusive  waves  described  in  the  Interim  Report  of  January  31, 
1983.  Further,  this  work  represents  a  major  advance  in  the  mathe¬ 
matical  structure  of  coupled  diffusion  equations. 

As  a  model  for  such  a  system  of  equations  we  consider  the 
transport  of  mass  and  energy  through  coupled  transport  equations. 
The  interaction  between  mass  and  heat  diffusion  (transport)  co¬ 
efficients  is  a  source  of  many  new  effects  described  by  the  follow¬ 
ing  prototype  system  of  equations : 

3tp  =  3xDl(p'T)3xp'*  (1) 

p3tT  *  3xpD2(p'T)3xT;.  %  <2) 


®1  B1  a2  B2 
D1  =  Dolp  T  '  D2  *  Do2p  T  ' 


where  3t  and  3X  are  partial  derivative  operators,  a,  8,  Dq1 

Dq2  are  constants  and  x  e[-l,l].  We  assume  slab  symmetry,  but  our 

results  are  valid  for  other  symmetries  as  well,  p  and  T  represent 
mass  and  temperature  respectively. 

The  initial  date  for  density  p(x,o)  and  temperature  T(x,o) 
are  defined  with  homogeneous  boundary  conditions  of  either  the 
first  type  (Dirichlet) , 

T(±l,t)  =  0  ,  p(±l,t)  =  0  ,  (4) 


or  the  third  type  (convective) , 

3xT  ±  hgT  =  0  ,  3xp  ±  hpp  =  0  at  x  =  ±1,  (5) 

where  h„  and  h  are  constants, 
e  p 

For  Dirichlet  boundary  conditions  [Eq.  (4)],  the  elliptic 
parts  of  Eqs.  (1)  and  (2)  are  degenerate  on  the  boundary,  and  a 
solution  may  exist  only  in  a  weak  sense.  This  singular  behavior 
of  the  equations  on  the  boundary  prevents  either  density  or  temp¬ 
erature  from  vanishing  when 

<  <*2  and  0^  <  $2  •  (6) 

When  this  occurs,  no  zero  Dirichlet  boundary  conditions  should 
be  iptposed  on  the  temperature. 

Similarly,  when 

*.*■ 

“l  >  °2  but  ^1  -  *2  +* l'  (?) 

the  density  must  be  bounded  but  cannot  be  required  to  vanish. 

When  *  a2  and  0^  <  02'  density  blocks  the  diffusion  of 'heat 
and  the  diffusing  plasma  becomes  isothermalized.  When  0^  =  02 

and  <»i  >  o2»  particles  cannot  escape  and  the  system  evolves  toward 
a  cold  constant-density  state. 

To  simplify  our  work  initially  we  have  concentrated  on  the 
convective  boundary  conditions  [Eq.  (5)]. 

From  previous  studies  of  a  single  diffusion  equation  des¬ 
cribing  diffusion  of  either  mass  or  energy,  we  expect  a  simple 
pattern  to  emerge  out  of  nearly  arbitrary  initial  data.  The 


separable  solutions  of  Eqs.  (1)  and  (2)  must  satisfy  special 
initial  conditions,  but  they  are  of  prime  importance.  Indeed, 
extensive  numerical  experimentation  has  shown  for  almost  arbi¬ 
trary  initial  conditions  that  the  system,  after  a  short  tran¬ 
sient  time,  either  evolves  toward  a  time-space  separable  solution 
or  actually  converges  to  it. 

To  summarize  briefly  our  results  so  far,  we  have  identified 
two  conceptually  different  diffusion  regimes.  In  the  first  re¬ 
gime  (a^  o2  and  B2  >  8^)  ,  every  initial  state  transits  into  a 
universal  diffusion  mode  given  by  the  space-time  separable 
solution.  The  decay  rate  of  this  asymptotic  solution  is  known 
a  priori  unles  a2^i  =  a1^2'  case  a  global  analysis  is 

needed.  In  the  second  regime,  diffusion  is  inhibited  and,  al- 

•  •»  *  % 

though  the  system  evolves  toward  the  separable  form,  it  cannot, 

mathematically  speaking,  attain  this  form.  The  nonlinear  inter¬ 
play  between  density  and  the  temperature  always  inhibits  the 
diffusion  of  either  temperature  or  density. 


3.  Exact  Integrals  of  the  Emden-Fowler  Equation 

The  celebrated  Emden-Fowler  equation  (henceforth  referred 
to  as  E.F.E.)  appears  in  various  branches  of  physics  and  en¬ 
gineering  and  as  such  was  and  still  is  a  subject  of  extensive 
analysis.  A  review  by  Wong  [1]  summarizes  the  investigations 
concerning  the  qualitative  properties  of  this  equation  and  its 
generalizations.  In  a  different  vein  Ames  and  Adams  [2]  employ 
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a  group  method  to  transform  the  E.F.E.  stated  as  a  boundary  value 
problem,  into  an  initial  value  problem  which  then  becomes  an 
easy  numerical  task. 

Our  interest  is  different;  we  focus  our  attention  on  an 


analytical  integration  of 

E.F.E.  written  as 

ii  •  i  i  rn  n 

xy"  -f  zy'  +  ax  y 

=  0 ,  a  =  const. , 

(la) 

or  in 

one  of  its  slightly 

generalized  forms 

xy"  +  (1+8)  y'  +  axmyn  =  0 

(lb) 

or 

<xv+Vr  +  axV 

=  0 

(lc) 

where 

m  =  1-a,  8  =  v-m. 

Related  to  our  interest  are  several  approaches  for  finding 

first  integrals  of  dynamical  systems  that  have  been  recently 

presented.  They  use  the  E.F.E. ,  (Eq.(la))  or  one  of  .its 

•  * 

variants  (Eqs.  (lb)  or  (lc))  as  a  test  case  to  demonstrate  the 
applicability  of  the  advocated  method.  These  procedures, 
generally  speaking,  are  either  variational,  group-variational 
(i.e.,  exploit  a  la  Noether  the  invariance  of  the  Lagrangian) 
or,  Ui_e  an  invariance  property  of  partial  differential  equations 
that  the  searched-after  first  integral  satisfies.  These  methods 
will  be  reviewed  and  compared  with  our  approach.  In  our  work, 
however,  we  approach  this  problem  directly  and  derive  two  suf¬ 
ficient  conditions,  of  which  at  least  one  appears  to  be  the 
hitherto  unknown.  When  either  of  these  conditions  is  satisfied, 
one  is  ensured  not  only  of  the  existence  of  first  integral  of 


■v  ..  *.* 
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motion  but  a  total  integration  of  E.F.E.  The  point  that  every 
first  integral  of  E.F.E.  may  be  brought  to  an  autonomous  form 
and  thus  further  integrated  is  trivial  but  was  surprisingly 
enough  unnoticed  in  previous  works.  It  is  noteworthy  that  each 
integrable  case  generates  a  one  parameter  family  of  integrable 
Emden-Fowler  equations. 


4 .  New  Results  in  the  Theory  of  Turbulent  Self-Diffusion 

The  work  described  in  this  and  the  following  section  is 
being  carried  out  by  E.  Hameiri. 

We  consider  the  diffusion  of  a  fluid  element  as  a  result  of 

stationary  and  homogeneous  turbulence.  Thus,  we  know  the  same- 

*  . 

time  statistics  of  the  turbulent  velocity  field  u  <u.  (0) -u. (0) > 

*»»  K  **  Xi 

2 

=  a^  and  would  like  to  determine  the  expected  deviation 

X (t)  of  a  fluid  particle  from  its  original  position. 

Y(t)  =  d  <|x(t)  -  x(0) |  >,  where  k  indicates  a  Fourier  mode  and 

*»* 

d  is  the  dimensionality  of  the  space. 

There  have  been  many  attempts  to  solve  this  problem.  For 
example,  Taylor  and  McNamara,  (Phys.  Fluids,  1973),  after  a 
number  of  approximations  obtained  the  result 


J?  Y  =  1  z  e-*/zk  ya2 

d?  d  k  k 


Y  (0 )  =  Y  (0 )  =  0 


which  determines  Y(t).  The  result  was  below,  but  rather  close 
to,  numerical  calculations.  This  work  however  had  some  short¬ 
comings.  First,  it  dealt  only  with  the  two-dimensional  case  d=2. 
Secondly,  it  required  some  precise  knowledge  of  the  flow,  namely 


that  the  vorticity  was  advected  by  the  fluid,  and  the  approxi¬ 
mation  making  use  of  this  information  appeared  only  one  of 
many  others  that  could  have  been  made. 


Our  approach  circumvents  the  need  for  knowledge  of  the 
flow  by  using  a  Lagrangian  variational  principle  which  determines 
it.  Using  then  a  well  known  (but  not  necessarily  valid)  ap¬ 
proximation  for  a  turbulent  medium  (based  on  non-dynamic  con¬ 
siderations)  as  an  additional  constraint,  we  get  the  result 


d2  22  l/2 

■A  Y  >  4(Za2)  (Ze 

dt^  ~  a  K 


-k2y,2, 

ak 


Inequality  holds  because  of  the  dropping  of  a  constraint  (the 

«  * 

knowledge  of  the  end  points  of -.all  particle  trajectories  in 
time) .  This  result  is  obtained  for  arbitrary  d,  and  is  in  strik¬ 
ing  contrast  to  Taylor's  relation.  A  Schwartz  inequality  .argu¬ 
ment  shows  our  Y (t)  to  be  larger  than  Taylor's.  Our  result  can 
be  improved  to  give  a  somewhat  closer  bound  by  using  additional 
information  on  the  flow  as  further  constraints  on  the  variation, 
e.g.,  conservation  of  enstrophy. 

This  novel  approach  to  the  diffusion  problem  can  be  similarly 
used  by  employing  other  variational  principles  which  determine 
the  same  flow.  In  particular,  one  can  get  a  lower  bound  on  the 
diffusion  by  using  a  principle  "dual"  to  the  first  one.  We  did 
not  succeed  yet  in  obtaining  a  lower  bound  apparently  close  to 
Taylor's  and  further  analytical  and  numerical  work  is  needed. 


5.  Mathematical  Theory  of  the  Essential  Spectrum  in  Maaneto- 


The  linearized  MHD  equations  present  a  rare  example  of 
differential  operator  which  has  a  non  pure  point  spectrum, 
even  though  it  may  be  defined  on  a  finite  domain  and  with 
coefficients  as  smooth  as  we  wish.  This  stands  in  marked  con¬ 
trast  to  the  much  discussed  Schrfldinger  operator.  From  a 
practical  point  of  view,  the  MHD  spectrum  determines  the  time 
evolution  of  small  perturbations  of  a  plasma  about  an  equili¬ 
brium  state.  The  singularities  in  the  spectrum  appear  in  the 
configurations  common  in  magnetic  confinement  experiments,  but 
the  underlying  causes  affect  the  behavior  of  all  magnetized 
plasmas,  e.g.,  solar  wind  plasma. 

The  present  work  represents  a  rigorous  mathematical  in¬ 
vestigation  of  the  subject.  Some  results  were  nevertheless 
known  previously,  usually  derived  by  heuristic  arguments.  In 
particular,  Grad,  Pao  and  others  derived  equations  which  de¬ 
termined  the  so-called  "Alfven"  system,  while  the  existence  of 
"ballooning  modes"  was  discovered  in  the  past  few  years.  These 
latter  modes,  which  were  derived  as  a  stability  criterion,  are 
shown  by  us  to  be  another  part  of  the  essential  spectrum.  (This 
spectrum  is  defined  as  the  whole  spectrum  except  for  all  discrete 
eigenvalues  of  finite  multiplicity.) 

The  derivation  of  the  essential  spectrum  is  done  by  us  by 
finding  an  "approximate  eigenfunction",  i.e.,  by  constructing 
a  sequence  of  functions,  not  converging  to  zero  in  the  norm,  such 


-13- 


that  the  eigenvalue  equation  is  satisfied  in  the  limit.  It 
always  happens  that  the  limiting  function  needs  to  be  localized 
in  space  near  a  single  field  line,  and  a  reduced  one-dimensional 
eigenvalue  equation  along  the  field  line  determines  points  in 
the  essential  spectrum  of  the  original  problem.  There  are  many 
possible  directions  of  strong  localization,  one  of  which  yields 
the  "Alfven"  spectrum  and  others  the  "ballooning  spectrum."  In 
the  case  of  axisymmetry,  the  spectrum  is  the  union  of  the  spectra 
corresponding  to  Fourier  modes  in  the  ignorable  direction,  and 
we  proved  that  each  such  part  consists  of  exactly  an  Alfven 
spectrum  plus  discrete  modes.  Thus  ballooning  modes  must  be 
the  accumulation  points  of  discrete  eigenvalues,  the  accumula¬ 
tion  occuring  over  non  empty  intervals.  Our  proof  here  utilized 
the  properties  of  compact  operators  in  the  theory  of  perturbations 

s  *  ^ 

»  . , 

of  the  spectrum  of  operators.  ' 

Finally,  we  gave  a  physical  interpretation  to  the  essential 
spectrum.  The  need  for  localized  eigenfunctions  indicates  the 
existence  of  one-dimensional  wave  propagation  along  magnetic 
field  lines.  Indeed,  one  such  wave  is  known  as  the  Alfven  wave, 
and  in  the  linearized  system  another  such  wave  exists.  Looked 
at  from  this  point  of  view,  solutions  of  the  localized  equations 
we  obtained  are  needed  in  order  to  determine  plasma  behavior 
in  all  configurations  with  a  magnetic  field,  e.g.,  whistler 
plasmas. 


AFOSR  -  Interim  Technical  Report 
January  1984 


List  of  Publications 


1)  Bayliss,  A.,  "An  Iterative  Method  for  the  Helmholtz  Equation," 

J.  Computational  Physics  £9  (3),  March  1983. 

2)  Bayliss,  A.  and  L.  Maestrello,  "On  the  Interaction  of  a  Sound 

Pulse  with  the  Shear  Layer  of  an  Axisymmetric  Jet  II. 

Heated  Jets,"  to  appear  in  Journal  of  Sound  and  Vibration. 

3)  Bayliss,  A.  and  L.  Maestrello,  "  Flow  and  Far  Field  Acoustic 

Amplification  Properties  of  Heated  and  Unheated  Jets,”  to 
appear  in  revised  form  in  AIAA  Journal. 

4)  Bayliss,  A.  and  E.  Turkel,  "Far  Field  Boundary  Conditions  for 

Compressible  Flows,  to  appear  in  J.  Comp.  Phys. 

5)  Bayliss,  A.,  E.  Turkel  and  C.  Goldstein,  "An  Iterative  Method 

for  the  Helmholtz  Equation,"  to  appear  in  J.  Computational 
Physics. 

6)  Grad,  H.,  "Correlations,  Fluctuations  and  Turbulence  in  a 

Rarefied  Gas,"  in  Long-Term  Prediction  in  Dynamics,  eds. 

C.W.  Horton,  Jr.,  L.E.  Reichl  and  A.G.  Szebehely,  John 
Wiley  &  Sons,  Inc.,.  1983.  *~ 

•,  . 

7)  Grad,  H.  and  S.L.  Wen,  "Equilibrium  and  Adiabatic  Compression 

of  a  Free  Boundary  Doublet,"  SIAM  J.  on  Appl.  Math.  £3  (1), 
February  1983. 

8)  Grossmann,  W.  and  J.  Teichmann,  "The  Effect  of  Stochastic 

Fluctuations  on  MHD  Waves,”  in  preparation  for  publication 
in  Phys.  Fluids. 

9)  Hameiri,  E.  and  H.A.  Rose,  "Magnetohydrodynamic  Fluctuations 

Near  Thermal  Equilibrium,"  Phys.  Fluids  £5  (12),  December  1982. 

10)  Hameiri,  E.,  "The  Equilibrium  and  Stability  of  Rotating  Plasmas," 

Phys.  Fluids  2£  (1),  January  1983. 

11)  Hameiri,  E.,  "Adiabatic  Compression  of  Rotating  Plasmas," 

Phys.  Review  A,  2 _7  (2),  February  198  3. 

12)  Hameiri,  E. ,  "Variational  Principles  and  Adiabatic  Compression 

of  Rotating  Plasmas,"  in  preparation. 


-2- 


13)  Hameiri,  E.  and  P.  Laurence,  "The  Ballooning  Spectrum  of 

Rotating  Plasmas,"  to  appear  in  J.  Math.  Phys. 

14)  Hameiri,  E.,  "On  the  Essential  Spectrum  of  Ideal  Magneto¬ 

hydrodynamics,"  in  preparation. 

15)  Imre,  K.  and  H.  Weitzner,  "Relativistic  Broadening  Near 

Cyclotron  Resonance,"  submitted  to  Phys.  Fluids. 

15)  Mond,  M.  and  H.  Weitzner,  "Stability  of  Helically  Symmetric 

Straight  Equilibria,  "Phys.  Fluids  25  (11),  November  1982. 

16)  Rosenau,  P.,  "A  Non-Linear  Thermal  Wave  in  a  Reacting  Medium," 

Physica  5D,  North  Holland  Publishing  Company,  1982. 

17)  Rosenau,  P.,  "Thermal  Waves  in  an  Absorbing  and  Convecting 

Medium,"  Physica  265  D,  North  Holland  Publishing  Company, 

1983. 

18)  Turkel,  E.  and  A.  Bayliss,  "Radiation  Boundary  Conditions 

for  Wave-Like  Equations,"  submitted  to  CPAM. 

19)  Weitzner,  H.,  "Linear  Wave  Propagation  in  Ideal  Magnetohdyro- 

dynamics,"  in  Handbook  of  Plasma  Physics,  eds.  M.N.  Rosenbluth, 
R.Z.  Sagdeev,  Vol.  I.,  Basic  Plasma  Physics  I,  edited  by 
A. A.  Galeev  and  R.N.  Sudan,  North  Holland  Publishing. Company, 

1983. 

20)  Weitzner,  H.,  "Geometry  Effects  on  Ion  Heating  at  Ion 

Cyclotron  Frequencies,"  Phys.  Fluids  26  (4),  '  ^rii  1983. 


